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ABSTRACT 

High spatial resolution observations show that high redshift galaxies are undergoing intensive evolution of 
dynamical structure and morphologies displayed by the Ha, H/3, [O III] and [N II] images. It has been shown 
that supernova explosion (SNexp) of young massive stars during star formation epoch, as kinetic feedback to 
host galaxies, can efficiently excite the turbulent viscosity. We incorporate the feedback into the dynamical 
equations through mass dropout and angular momentum transportation driven by the SNexp-excited turbulent 
viscosity. The empirical Kennicutt-Schmidt law is used for star formation rates. We numerically solve the 
equations and show that there can be intensive evolution of structure of the gaseous disk. Secular evolution of 
the disk shows interesting characteristics that are 1) high viscosity excited by SNexp can efficiently transport 
the gas from lOkpc to ~ lkpc forming a stellar disk whereas a stellar ring forms for the case with low viscosity; 
2) starbursts trigger SMBH activity with a lag ~ 10 8 yr depending on star formation rates, prompting the joint 
evolution of SMBHs and bulges; 3) the velocity dispersion is as high as ~ 100 km s" 1 in the gaseous disk. 
These results are likely to vary with the initial mass function (IMF) that the SNexp rates rely on. 

Given the initial mass function, we use the GALAXEV code to compute the spectral evolution of stellar 
populations based on the dynamical structure. In order to compare the present models with the observed 
dynamical structure and images, we use the incident continuum from the simple stellar synthesis and CLOUDY 
to calculate emission line ratios of Ha, H/3, [O III] and [N II], and Ha brightness of gas photoionized by young 
massive stars formed on the disks. The models can produce the main features of emission from star forming 
galaxies. We apply the present model to two galaxies, BX 389 and BX 482 observed in SINS high-z sample, 
which are bulge and disk-dominated, respectively. Two successive rings independently evolving are able to 
reproduce the main dynamical and emission properties of the two galaxies, such as, BPT diagram, relation 
between line ratios and Ha brightness. The observed relation between turbulent velocity and the Ha brightness 
can be explained by the present model. High viscosity excited by SNexp is able to efficiently transport the gas 
into a bulge to maintain high star formation rates, or, to form a stellar ring close enough to the bulge so that it 
immigrates into the bulge of its host galaxy. This leads to a fast growing bulge. Implications and future work 
of the present models have been extensively discussed for galaxy formation in high-z universe. 

Subject headings: galaxies: evolution - galaxies: high-redshift 



1. INTRODUCTION 

Star forming galaxies at redshift z ~ 1.5-3.5 around the 
peak of cosmic star formation and quasar activity have a large 
variety of dynamical properties and morphologies (Forster 
Schreiber et al. 2006; Genzel et al. 2006, 2008; Wright et 
al. 2007, 2009; Bournaud et al. 2008). Four main char- 
acteristics different from the local galaxies are found from 
the large available high-z samples (Genzel et al. 2008; 
Forster Schreiber et al. 2009; Lehnert et al. 2009; Lemoine- 
Busserolle & Lamareille 2009): 1) higher velocity dispersion 
Vrot/c < 10 throughout entire regions of galaxies at high red- 
shift than V rot /a ~ 15-20 in local galaxies (Dib et al. 2006), 
where V 10t is the rotational velocity and a is the velocity dis- 
persion; 2) high-z galaxies are more irregular and asymmetric 
in shape, which is not related with ongoing mergers (Shapiro 
et al. 2008); 3) they show highly clumpy disks, typically with 
a scale of < lkpc (e.g. Genzel et al. 2008); 4) very high sur- 
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face brightness of the recombination lines (e.g. Lehnert et al. 
2009). There is no doubt that high-z galaxies are undergoing 
most violent formation and evolution. 

The "heating" mechanism of high velocity dispersion of 
clumps remains open, which is thought potentially to be re- 
sponsible for the transportation of angular momentum of gas. 
Turbulence driven by self-gravity instability is not favored by 
evidence that such flows are not able to produce the dispersion 
and emission via mergers (Lehnert et al. 2009). However, 
there is increasing evidence for the turbulence excited by su- 
pernova explosion (SNexp) of young massive stars during star 
formation from both numerical simulations or theoretical ar- 
guments (Wada & Norman 2001; Mac Law & Klessen 2004; 
Wang et al. 2009b) and observations (Dib et al. 2006; Lehn- 
ert et al. 2009). The strong correlation between the Eddington 
ratios and specific star formation rates in type 2 AGNs implies 
the roles of SNexp in exciting turbulence, which triggers an 
accretion flow to central SMBHs (Chen et al. 2009; see also 
Watabe et al. 2008). In particularly, the strong correlation be- 
tween the velocity dispersion and star formation rate density 
in both local (Dib et al. 2006) and high-z galaxies (Lehnert et 
al. 2009) indicates the role of the star formation in triggering 
turbulence. We have to note also that turbulence can enhance 
star formation in term of the compressed density by shocks 
(Silk 2005). Anyway, SNexp inputs energy into interstellar 
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medium causing turbulence somehow. 

It has been realized that the star formation itself plays a 
key role in the formation of galactic structures. The well- 
known exponential disks of most S and SO galaxies (e.g. Free- 
man 1970; Wevers et al. 1986) can be explained by that the 
gaseous disk is evolving like a viscose accretion disk, and 
the exponential structure reaches when the infalling timescale 
equals to the star formation (Lin & Pringle 1987; Lu & Cheng 
1991), or by self-regulations of star forming galaxies (Dopita 
& Ryder 1994; Silk 1997). Tracing back the galaxies to high 
redshift, they must have been undergoing intensive star forma- 
tion and similar to those galaxies at high redshift which have 
exponential disks today. Wang et al. (2009b; hereafter Paper 
I) set up equations describing the self-regulating gaseous disk 
driven by SNexp feedback for high redshift galaxies, and, for- 
tunately, find that such a set of complex equations has analyt- 
ical solutions if the star formation law is a linear relation with 
gas density. Paper I shows the important roles of SNexp in 
the secular evolution to supply gas to ~ lkpc regions and find 
that the prominent feature of star forming galaxies in high-z 
universe is a natural consequence of turbulent viscosity driven 
by SNexp. Very recently, Kumar & Johnson (2010) extend the 
study of the roles of SNexp in the transportation of fueling gas 
at ~ 1 kpc scale and confirm that SNexp is an efficient way to 
transport angular momentum outward. The advantage of this 
self-regulation is that no external torque is needed to transport 
gas into inner region. It is obvious that the dynamical struc- 
ture and morphology is highly time-dependent, which should 
be given by a self-consistent way. 

On the other hand, coevolution of SMBHs and galaxies has 
been generally formulated by the correlations between SMBH 
mass and bulge luminosity (Magorrian et al. 1998); or dis- 
persion velocity (Ferrarese & Merritt 2000; Gebhardt et al. 
2000). Great attempts of phenomenological models have been 
made (e.g. Granato et al. 2001; Colpi et al. 2006), however, 
dynamics of the coevolution is poorly understood so far. On 
the side of theoretical expectation, gaseous clumps in galaxies 
not only undergo star formation, but also will be partially de- 
livered into the central regions of galaxies by SNexp (Paper I). 
As a natural consequence of evolution of clumps, bulges are 
growing and SMBHs do subsequently. The time lag of SMBH 
activities relative to the starbursts is determined by the star 
formation itself. In principle, solutions of the time-dependent 
equations describing the dynamical structures of gaseous disk 
can present the whole story of the coevolution. 

In this paper, we extend studies of Paper I in detail for the 
intensive evolution of dynamical structure and morphologies 
of high-z galaxies. In §2, we use the Kennicutt-Schimidt 
law to calculate the star formation rates, and incorporate the 
rates into the dynamics of the gaseous disk through the SNexp 
excited-turbulence viscosity. Growth of bulges and SMBHs 
are discussed as well as the time lag of triggering SMBH ac- 
tivities with respect to starbursts. Photoionization model is 
presented in §3 for different geometries and conditions in the 
evolving gaseous disk. Detailed applications to two high-z 
galaxies are given in §4. We extensively discuss implications 
and future improvement of the present model in §5. Conclu- 
sions are drawn in the last section. 

2. DYNAMICS AND EVOLUTION OF THE GASEOUS DISKS 

2.1. Dynamical equations 

Spatial resolution of Integral Field Spectroscopy is ~ 
O.lkpc at z ~ 2, which is the maximum size of clumps re- 



solved observationally (Genzel et al. 2008; Jones et al. 2010). 
This size is still smaller than the characterized length of galax- 
ies (~ 10 kpc), but the number density of the clumps are dense 
enough to approximate the system of clumps as continuum 
fluid. It is valid to use the equations of continuum medium, 
but we treat the gas as clumps in the photoionization model. In 
Paper I, we derive a series of equations to describe the secular 
evolution of a gaseous disk by including feedback of star for- 
mation. We assume that turbulence excited by SNexp of mas- 
sive stars is driving the angular momentum outward, and mass 
dropout is happening due to ongoing star formation. Actually 
the role of supernova-driven turbulence has been stressed by 
evidence of correlation between the dispersion velocity and 
star formation rates (Dib et al. 2006; Lehnert et al. 2009). 
For a convenience, we list the equations here. With the mass 
dropout, the mass conservation equation reads 



R- 



<9E a d 



+ —{RV R Z g )+R^ = 0, 



dt dR 



(1) 



where Vr is the radial velocity of the gas, E g is the surface 
density of gas, and E* is the surface density of star forma- 
tion rates. All the parameters E, E*, Vr, f2, Q, P lm h, Vturb, H, 
and p g are functions of radius R unless explanations. We ne- 
glect mass injection of stellar winds back to the gaseous disk 
or escapes of winds from galaxies. Star formation process re- 
moves some angular momentum of gas, thus the conservation 
equation of angular momentum is given by 

Rj t (x s R 2 n) + ^ (Rv R z g R 2 n) - ±- 1| +R 3 n£+ = o, (2) 



dR 



2tt dR 



where Q = -27r/?VE g (dVl/dR) is the viscosity torque, v is 
kinematic viscosity, and £1 is the angular velocity. Combining 
equations (1) and (2), we have 



<9E e 



1 d 



dt 2-kR dR 



d (R 2 n) 

dR 




(3) 



Let X = R 2 Vt, equation (3) is converted into a more concise 
form 

dt InRdRdX 1 - g ' K> 

where the Kennicutt-Schmidt (KS) law as E* = AEg is used 
and A is a constant (Kennicutt 1998fl This equation includes 
the star formation process which has two major effects: con- 
sume the gas and apply energy to make turbulence by SNexp. 
Imposing adequate boundary and initial conditions to equa- 
tion (4), it will produce the structure of the gaseous disks. 

To further simplify equation (4), we have to specify the 
relation between E g and Q through the viscosity parame- 
ter v. Thermal turbulence (currently thought as origin from 
magnetohydrodynamic -rotation instability in Balbus & Haw- 
ley 1998) is assumed to make viscosity for transportation of 
gas angular momentum as known as the a-prescription in 
standard disk model (Shakura & Sunyeav 1973). Obviously, 
the clouds in the star forming regions are not hot enough to 
trigger thermal turbulence. However, the turbulent velocity of 

5 Star formation law at high— z could be different from the KS law (Gnedin 
& Kravtsov, 2010), but we ignore the difference for simplicity. The influence 
of the star formation law on the structure of the gaseous disk can be generally 
estimated by equations (15-17) and (19-21). See a brief discussion in §2.2.6. 
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TABLE 1 

Summary of all input parameters in the model 



Parameter 


Units 


Physical meanings and values 


C* 


Gyr 1 


the star formation efficiency (~ 0.25) 


^SN 


ergs 


the kinetic energy of supernova explosion (~ 10 51 ) 


fsN 


M ~A 


the parameter of converting SF rates to SNexp rates 


V c 


km s -1 


a characterized velocity of the rotation 


R, 


kpc 


the character radius with velocity V c 


Ro 


kpc 


the initial radius of the gaseous ring 


So 


M Q pc" 2 


initial surface density of the ring 


M 


1O 1O M 


mass of the initial gaseous ring 


a 




the viscosity parameter 


7 




the index of the KS law (~ 1 .4) 


? 


M~ l 


efficiency of the supernova explosion 


A R 


kpc 


the width of the initial ring 


At 


Gyr 


the interval of mergers (~ 0. 1) 


f 




a ratio of radiation to turbulent pressure in photoionization 


«H 


cm -3 


density of ionized gas 



^100 km s -1 much higher than thermal sound speed can effi- 
ciently make the turbulent viscosity. This has been supported 
by that turbulence excited by SNexp is strong to power dy- 
namical viscosity and transport the angular momentum (Wada 
& Norman 2001; Chen et al. 2009; Paper I and Lehnert et al. 
2009). We modify the form of a- viscosity prescription via 
replacing the sound speed by the turbulent velocity driven by 
SNexp 

v = aV tm H, (5) 

where a is a constant, V tur is the turbulence velocity driven 
by SNexp and H is the thickness of the disk. This prescrip- 
tion can be explained by that clouds are communicating with 
each other through the turbulent velocity, and the maximum 
scale of turbulence can not exceed the height of the disk. All 
the uncertainties of the turbulent viscosity are absorbed into 
the parameter a. It is thus expected a < 1 in the turbulent 
disk. We take a as a constant throughout the paper. Detailed 
transportation of angular momentum driven by SNexp can be 
found in Rozyczka et al. (1995) and briefly in Collin & Zahn 
(2008), but we use the simplified version. 

Assuming the kinetic energy is channeled into turbulence, 
energy equation of turbulence excited by SNexp is given by 



Pg^tur _ PgVt 



g ' tur 

H 



■■ eS+Esn, 



(6) 



where f^is = H /V tm: is the dissipation timescale of the tur- 
bulence, e is the efficiency converting kinetic energy of 
SNexp into turbulence, 5* is the SNexp rate and £sn is the 
SNexp energy. SNexp rate strongly depends on the ini- 
tial mass functions [IMF; N(M*)], we have a factor /sn = 
Jm™* N(M*)dM*/M tot to produce SNexp in a time, where 

M tot = f^ m "N(M*)dM*. We then have 5* = fsx£*/H. It is 
convenient to introduce £ = e /sn for discussion (see §2.2.5 for 
details). Considering relations of E g = p g H and the KS law, 
we have 

V 3 

* Hi 



(7) 



We note that this equation energizing the turbulence is highly 
simplified, especially the efficiency £. The dependence of £ on 
the IMF is discussed in §2.2.5 and §5. Additionally, SNexp 
usually lasts ~ a few 10 4 years (depending on its surround- 
ings) divided into 3 phases, which is much shorter than the 



typical timescale of dynamical evolution (<~ 10 7-8 yrs). We 
thus neglect the time-dependent effects of the SNexp. 
We assume an equilibrium in vertical direction 



gravj 



(8) 



where P tm = p g V 2 m is the turbulence pressure and P gmv is the 
gravity in vertical direction. There are two cases to balance 
the turbulence pressure (Binnery & Tremaine 2008; Vollmer 
& Beckert 2002): A) local stars provide the vertical force 
and B) central mass dominates. The two cases have differ- 
ent vertical equilibrium equations. In Case A, p g <C and 
M g (R) <C M+(R), where p g (or p*) is gas (or stellar) density 
and M g (R) [or M+(R)] is enclosed mass of gas (or stars) within 
a radius R, namely, surface density of existed stars dominates. 
In such a case, we have Pg TaY — TrGS^Sg and V£ t = nGEJt 
(Binney & Tremaine 2008, Vollmer & Beckert 2002), 



Vro, 



H 

R' 



(9) 



where V rot is the rotation velocity at radius R and V tur is the 
turbulent velocity. Such a relation has been suggested by 
Genzel et al. (2008) to describe a large thin disk for their 
SINS galaxy sample. In Case B, M g (R) < 0.5(H /R)M(R), 
the disk is dominated by the central mass within R. Then 
it holds Pgrav = p g n 2 H 2 (Pringle 1981, Vollmer & Beckert 
2002). Combining V? ol = il 2 R 2 and S g = p g H, equation (8) 
is rewritten by 



V m t 



H 
R' 



(10) 



which is used to describe a thick or compact disk (Genzel et 
al. 2008). The different equilibrium in vertical direction will 
lead to different structures of the gaseous disk globally. It is 
worth pointing out that the vertical equilibrium is sensitive to 
£ as shown in subsequent sections. 

We specify rotation curves of the gaseous disk in order to 
solve equation (4). Following Binney & Tremaine (2008), we 
assume that galaxies have a logarithmic potential leading to 
the following rotation curve 



Vro,= 



l + (Rc/R) 



(11) 
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where R c is a critical radius and corresponds to the velocity 
y rot = V c /V2, and V c is the velocity at infinity. The rotation 
velocity and its gradient are then given by 



n = 



R 



V r R 



VR 2 +R 2 c' dR (r 2 +r2) 



3/2- 



(12) 



We would like to point out that the rotation curves are very 
similar for the mid-plane gaseous for both the disk-dominated 
and bulge-dominated. R c is a critical radius within which 
bulge's potential dominates, and it thus approximates to the 
size of its bulge. The present model focuses on the dynamics 
of the gas at the mid-plane, we thus adopt the same rotation 
curves for the two cases for simplicity. 
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FIG. 1. — Rotation curves for different V c and R c . Such a rotation curve is 
very similar to the observed in high-z galaxies (Genzel et al. 2008). Here we 
plot the curves for BX 389 and BX 482. The dotted line is the rotation curve 
used in the calculations for a general case, where we use V c = 300km s" 1 and 
R c = 3kpc. 

We list all the equations (4, 5, 7, 8, 9, 10) describing the 
evolution and structures of the gaseous disk with feedback of 
star formation. We stress that the above equations are actually 
averaged in vertical direction and the parameters obtained in 
these equations are the values at the mid-plane of the disk. 
Calculations of vertical structures of the disk are definitely 
important to compare with observations, but it is beyond the 
scope of the present paper. 

Initial and boundary conditions should be imposed to equa- 
tion (4) for the disk. The simple initial condition is a pulsar 
injection with a Gaussian profik^J 



S„(i?,0): 



TinR, 



0"0 



2ttA s 



exp 



(R-Ro) 2 



2A 2 



(13) 



at radius Rq and with a width A s as a consequence of a mi- 
nor merger (Hernquist & Mihos 1995), or cold flow (Dekel 
et al. 2009). Though it is very simple, numerical simulations 
of Hernquist & Mihos (1995) and Dekel et al. (2009) show 
the reality of the initial condition (equation 13), which is also 
supported by high-z galaxies (Genzel et al. 2008). Integrat- 
ing the ring, we have its total mass Mq = J 2nY, g RdR, which 
is the merged gas mass from the merged galaxy. The Gaus- 

6 Actually, it allows us to take any form of the initial profile of the gaseous 
ring in the numerical calculations. The Gaussian profile as initial condition 
is oversimplified, however, this form makes evolutionary properties of the 
gaseous ring more clear. For more extended initial gaseous rings, their evo- 
lution is approximately composed of multiple Gaussian rings. The behaviors 
can be found in Wang et al. (2009). 



sian profile used here is not physically necessary, but only 
for simplification to produce the main features of the gaseous 
disk. Any form of initial conditions can be performed by the 
numerical scheme in the Appendix. 
The inner and outer boundary conditions are given by 



£(0,f) = 0, and 



dG(X,t) 



dX 



= M(X max ) = 0, (14) 



x=x„. 



where X max is the maximum of the parameter X. The inner 
boundary is given by the torque-free condition whereas there 
is no steady injection to galaxies at the outer boundary. 

The viscosity parameter (a) is a free one, but it is still 
constrained by the simple argument below. The character- 
ized length of medium influenced by SNexp can be estimated 
by the fact that the SNexp stops sweeping the medium un- 
til its kinetic energy is exhausted by the medium. Since 
the kinetic energy of SNexp is equal to the thermal energy 
of swept medium, we have £sn ~ A%R\„nm p /3, namely, 

Rm ~ 0.2 (£51 / 'njTj) 1 3 kpc, where n-i = «/10 2 cm -3 is the 
number density of the medium, T2 = T/10 2 K is the temper- 
ature and £51 = £sN/10 51 erg. This length is comparable with 
the turbulence size due to SNexp though the actual situation is 
very complicated (e.g. Vollmer & Beckert 2002). For a thick 
disk with height of ~ lkpc, we have a ~ Rsn/H ~ 0.2, which 
absorbs the uncertainties of viscosity, dissipation length and 
turbulence driven by SNexp. We use a = 0.2 throughout the 
present paper. This is also important for calculations of pho- 
toionization, where we find /?sn is naturally about the size of 
clumps (see equation 32). 

In this paper, we use the empirical KS law = AS? , where 
7= 1.4, A = 2.5 x 10" 4 , and E g are in units of M yr _1 kpc~ 2 
and M Q pc~ 2 , respectively. For simplicity in making equation 
(4) dimensionless, using = (10~ 5 A)yr _1 , we rewrite the KS 
law as = CtSgSJl , where S gi i = E g /M Q pc" 2 , for the nu- 
merical scheme. 

In a summary, we list all the necessary equations and con- 
ditions in this subsection. All the parameters invoked in the 
models are given in Table 1 . We fix several parameters, such 
as C*, £sn, a > At and 7 throughout the paper. We do not 
study the cosmic evolution of galaxy in the present paper, but 
focus on the evolution at dynamical timescale. The rotation 
curves are then fixed in the evolution. We stress the differ- 
ences of the present equations from that in Paper I, in which 
linear star formation law and constant rotation curves are used 
for analytical solutions. Only numerical solutions exit for the 
present cases. We describe the numerical scheme to solve the 
equations in the Appendix. 

2.1.1. Case A: Stellar mass dominated disk 

For a fast glance of dependence of the solution on the pa- 
rameters, we have form solutions of the gaseous disk from 
combining equations (7) with (9) for rough estimates of the 
structure of the gaseous disk before numerically solving the 
equations. They are 



y tur =ii3.8e-2£5iC_ 10 s: 



7-1 



tfioV 30 2 okms \ 



27-2 r>3 



W kpc, 



v = 33.4 «o.2^ 2 £ 5 3 i^ 10 £2 7 " 3 <)^oo kpc 2 Gyr" 1 , 



(15) 



(16) 



(17) 
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FIG. 2. — Structures and evolution of gaseous disks for Case A. From the upper to bottom panels are gas surface density, star formation rate surface density, 
stellar surface density, height, turbulent velocity and the Toomre parameter, respectively. This plot is for the given initial conditions with So = 50Mqpc , 
(corresponding to S* o = O.O6M0kpc~ 2 yr~'), Rq = lOkpc, A.r = 0.5kpc, Mq = 3.25 X 1O 1O M0. Lines with different colors (from red, green, blue, to pink) 
refer to different times from 10 6 yr to 10'yr with an interval A logf = 0.25. And the black one is refer to the initial condition. 



R/10 kpc, E 51 = £ S N/10 51 erg, £_ 2 = 
V300 = Vrot/SOOkms" 1 , 



: C*/2.5 x 1(T 10 yr" 1 



where ao.2 = a/0.2, Rio : 
e/10- 2 A^,C_i 0: 
E 2 = S g /2OOM pc- 2 . 

We find from the form solutions that: 1) the velocity dis- 
persion excited by SNexp can be comparable with the obser- 
vations; 2) the disk is thin as H/R ~ 0.1; 3) the infalling (or 
advection) timescale of gas is of f a d v ~ R 2 jv = R 2 w /v^ ^3.0 
Gyr due to the viscosity, where 1/33 = ^/33kpc 2 Gyr" 1 ; 4) V lul , 
H and v are very sensitive to the rotation velocity V rot , which 
fully depends on the host galaxies. Clearly the evolution of 
gaseous disk is sensitive to the host galaxies. 
With the rotation curve given by equation (11), we get v = 



^37-3 



R 2 + R 2 ) 



/R 4 and g = 27rC A V c E^- 3 E g (^ 2 + « 2 ) 5 / 2 , 



where ( A = a(£EsNC*) 3 /K- 8 - We make the equations dimen- 
sionless through r = fC*, r = R/R*, r c = R c /R+, x = X /(V C R^), 

= 1, and F A = 



where R* = R 2 max / \/R 2 + R 2 mlix making x n 
g/(2n( A RlV c M G pc- 2 ) = ^- 2 (r 2 + r 2 ) 5 / 2 . 
tion (4) as 



We rewrite equa- 



<9 2 F A A dY A , 



(18) 



where the coefficients are 



kf(x,Y A ) = 



and 



k£(x,Y A ) 



C,(r 2 + r 2 )^ ^ 
C4-R 2 (37-2)(2r 2 + r 2 ) A 

C(r 2 + r 2 )^ A 



Q A Rl(2rl + ri) A 



The initial condition (equation 13) can be converted 

into Y A (x,0) = (r 2 + r 2 ) 5/2 E^" 2 (r,0), where £ g ,i(r,0) = 

£o,in)/v / 27rcrexp [(r-ro) 2 /2<7 2 ] . The boundary condi- 
tions (equation 14) are converted into F A (0,f) = and 
dY A / dx{x maX: t) = 0. Equation (18) is a non-linear 2nd par- 
tial differential equation. Numerical scheme for the equation 
is given in the Appendix. 

2.1.2. Case B: Central mass dominated disk 
Combining equations (7) and (10), we have 

V tur = 184.7 £iEg£ 10 



0.5p0.5 r 0.5 yO.57-0.5 ff 0.5 v -0.5 Vrn -1 
-2 £ 51 ^-10^2 "10 v 300 Km s ) 

H = 6.2 eM^^MiV^ 5 kpc, 



(19) 
(20) 
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FIG. 3. — Structures and evolution of gaseous disk for Case B. We use the same values of the initial conditions as in Fig. 2 for case A, but the parameter £ is 
different. The global properties of the gaseous disk is similar for Case A and B, but they are more sensitive to £ for Case B. 



v = 232.5 a . 2 UE5iC- W Y? 2 - l R 2 m V^ Q kpc 2 Gyr" 1 . 



(21) 

Obviously, the structures of the Case B disks are quite dif- 
ferent from that of Case A: 1) The disks have higher velocity 
dispersion as high as 100 km s" 1 ; 2) a quite thick geometry 
with H/R ~ 0.5; 3) the infalling timescale is much shorter 
than that of Case A by a factor of ~ 7. Moreover, unlike Case 
A, the disks are not so sensitive to the rotation curves, how- 
ever, it still depends on the rotation curves. 

We use the similar parameters to make equation (4) dimen- 
sionless. From the rotation curve, we get v = Cfi^gT 0^? + 
R 2 ) and Q = 27rC B V c 2^~ 1 £ g fi 4 (fi 2 + R 2 T l/2 , where ( B = 
a^E SN CjV 2 . Using Y B = Q /(2^Cb^V c M oPC - 2 ) = r\r 2 + 
^'^Jj.we get 



d 2 Y } 



(22) 



where 



and 



*?(jc,y B ) = 



C,(r 2 c +r 2 ) — 

4 

Cb7(2a-2 + r 2 )rn 



2x2 



*?Cx,y B ) = 



CAr 2 + r z ) 
0,(2r? + rV J 



The initial condition is F B (x,0) = E^(r,0)r 4 (r 2 + r 2 ) -1 / 2 , and 
the boundary conditions are identical with Case A. 

2.2. Structure and evolution of a single gaseous ring 

We tested the present numerical scheme in Appendix by 
setting 7 = 1 to check if the numerical solutions are identi- 
cal to the analytical presented in Paper I. It turns out that the 
numerical scheme in the Appendix works very well for the 
non-linear differential equation. 

Figure 1 shows several different rotation curves determined 
by the two parameters V c and R c . These curves are character- 
ized by two parts: 1) the inner part, which has a linear increase 
with radius; and 2) the outer, which becomes quite fiat. The 
structure of the outer part (beyond R c ) is quite similar to the 
analytical solutions, however, in the inner part most of the gas 
will be exhausted for star formation in the gaseous disks with 
high and intermediate viscosity. This region is very important 
to supply gas to the central supermassive black holes. 

We set typical values of So, Rq, Rc an d V c to solve equations 
(18) and (22), respectively. Though the initial values of the 
parameters influence the properties of the gaseous disk some- 
how, the main driver of the solution is obviously the param- 
eter £, which represents the feedback strength of star forma- 
tion. Additionally, KS law is a non-linear relation and leads 
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to some interesting dependence of the behavior of the disk on 
the initial surface density. We focus on the ^-dependence of 
the disks for Case A and B and discuss the non-linear effects 
of the KS law. Table 2 gives values of £ used in calculations. 

2.2.1. stellar and gas surface densities 

Generally, there are three types of solutions for Case A and 
B found in the £ g -panels of Figure 2 and 3. For a disk with 
low-viscosity, as shown by Figure 2c and 3c, the gas ring 
slowly diffuses spatially only. This behavior is very similar to 
the classical gaseous ring described by Lynden-Bell & Pringle 
(1974). It forms stars around the initial radius appearing as a 
stellar ring or a stellar disk in the diffuse timescale, which is 
longer than that of the star formation. Infalling gas advected 
by the turbulence-viscosity driven by SNexp is not powerful 
enough to transport gas inward. Here we simply ascribe the 
low viscosity to the low SNexp rates, which could be caused 
by a very steep IMF. In such a case, the fraction of young mas- 
sive stars are fewer than that of flat one, which determines the 
SNexp rates and hence £. We also note that the gas is not able 
to diffuse into R c . The gaseous disk is mainly dynamically 
controlled by the part with constant rotation. The solution of 
the low-viscosity ring is very similar to the analytical one in 
Paper I. 




0.5 1 0.5 1 



R/R R/R 

FIG. 4. — The ratios of advection to star formation timescales at different 
radius and time, q < 1 means that advection of gas by the SNexp-driven 
viscosity dominates over the star formation. The panel a, b and c correspond 
to that in Figure 2 and 3. 

For a gaseous ring with high-viscosity as shown in Fig- 
ure 2a and 3a, the gas is efficiently delivered inward by the 
turbulence-viscosity torque through SNexp, forming a peak 
and a valley of the envelope of the £ g -evolution track at /? pea k 
and /? V aiiey> respectively. The ring has three phases during 
its evolution. First, it dramatically spreads into a thick disk 
within about a few 10 6 yrs. Second the ring keeps a thick disk 
within a few 10 7 yrs. And finally it goes into the last phase, in 
which the disk evolves into a stellar disk within a few 10 8 yr. 
Contrary to a ring with low viscosity, the high viscosity ring 



forms a stellar ring or a stellar disk far away from the initial 
radius of the gaseous ring by a distance of ~ 0.1 initial radius. 
The high-viscosity arises from higher SNexp rates, which im- 
ply a high star formation rate or a top-heavy IMF. 

For a gaseous disk with intermediate-viscosity as shown in 
Figure 2b and 3b, star formation is competing with transporta- 
tion of gas inward, showing a quite flat distribution of gas den- 
sity without peak and valley of the envelope of £ g -evolution 
track. We note that the properties of the £ g -evolution track 
depend on the rotation curves, the initial gas density and ra- 
dius, but we focus on understanding the role of star formation 
in S g -evolution track. 

The peak and valley of the S g -evolution track are de- 
termined by the roots of the equations of dH g /dR = and 
<9£ g /<9f = 0. We have R pel ± and R va n ey for the analytical so- 
lutions in Paper I. The present properties of S g -evolution 
tracks are similar that in Paper I, but the gaseous disk shows 
faster evolution than the analytical solution. This is be- 
cause the empirical star formation KS law is steeper than 
the linear relation. Considering the timescale of star forma- 
tion t* = £ g /£* = £g"7/C* = 06 S2 0,4 Gyr, and gas advection 
timescale f adv =R/V R ^ R 2 /v = 1 / C^R 2 Y? g J~ l) for case A and 

?adv = f° r case B> the solutions can be qualitatively 

understood through the ratio 

( 6.2 aoiCioCl-Ssi^To^oo 5 ^ 27 > 
q(R,t) = -f^=\ (23) 

* { 0-9 aJ 2 &EjVi 00 , 

for Case A and B, respectively. It is interesting to find that q 
just follows y rot for Case B. 

Figure 4 shows the parameter q(R,t). The sharp peaks in 
these panels correspond to the point (R zao ) where the drift 
velocity (Vr) is very low and close to zero. Beyond ^ zero , 
the gas is transported outward whereas within R zero the gas is 
advected inward. We would like to point out that only tiny 
fraction of the gas is transported outward (beyond the initial 
radius Rq) as shown in the S g -panels. It is clear to show from 
o-panels that the low-viscosity leads to q > 1, namely, star 
formation dominates whereas the high-viscosity does q < 1 
with dominant advection of gas. 

2.2.2. geometry and turbulence velocity 

We plot the height of the disks and the turbulence velocity 
in H— and V tur -panel in Figure 2 and 3, respectively. Gener- 
ally, H/R< 1 holds for Case A and B, and the disk at the final 
stage has H/R ~ 1. The gaseous disk becomes very sharp at 
the inner edge for the cases with high and intermediate viscos- 
ity. The sharp peak of the disk height at its inner edge arises 
from that the left gas is bloated by the intensive SNexp during 
the star formation in light of a relatively high £. For a ring 
with low viscosity, it simply spreads over the space with a 
constant H/R for all the time. This is very similar to the clas- 
sical solution of diffusion problem for a ring with a constant 
viscosity (e.g. Lynden-Bell & Pringle 1974). 

We plot the turbulence velocity in Figure 2 and 3. We find 
that V tm holds a rough constant over the disk, but evolving into 
smaller one with time. This is caused by the decreases of star 
formation rates. The low viscosity arisen by low SNexp rates 
only supports a thinner disk. This predicts a relation between 
the height and dispersion velocity, however, it is very difficult 
to test the relation since the height of the gaseous disks is 
unknown from observations. 
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2.2.3. Self-gravity 

Star formation happens in the self-gravity-dominated re- 
gions. The numerical solutions allow us to test the self-gravity 
of the gaseous disk through the Toomre parameter defined as 

2=-^r (24) 

where n = 40 2 + dfl 2 /d InR, c s is the sound speed, we plot 
Q in g~P ane l of Figure 2 and 3. The parameter Q actu- 
ally represents the ratio between the vertical gravity to the 
self-gravity. Q < 1 implies a self-gravity dominated solution 
whereas Q > 1 means the self-gravity can be neglected. We 
find Q is smaller than unity in most regions during the evo- 
lution, making the star formation ongoing. This is in agree- 
ment with observations of Genzel et al. (2008) showing evi- 
dence for globally unstable of disk. We would like to stress 
that the evolution of gaseous disk is controlled by the SNexp 
through viscosity, rather than by adjusting the Toomre param- 
eter Q <~ 1. This is different from that in Thompson et al. 
(2005) who suggest that Q ~ 1 keeps by adjusting the star 
formation rates. However, the disagreement between Wada 
& Norman (2001) and Thompson et al. (2005) could be re- 
moved if the gaseous disk is clumpy, in which the definition 
of the parameter Q should be modified below. 

We note that Q' = V tur K/nGY, g is widely used for clumpy 
disk in Vollmer & Beckert (2002) and regarded as a free pa- 
rameter. Actually Q' is very different from Q. Communica- 
tion between clumps or clouds is in term of the turbulent speed 
V tur , but the gravitation instability develops through location 
sound speed c s . Q' describes the relative strength of the tur- 
bulent ram pressure to the vertical direction self-gravity, and 
does not deliver the classical meanings of the Toomre param- 
eter. Unlike Thompson et al. (2005), given the star forma- 
tion law, Q' can be derived, which plays an essential role in 
clumpy disks. 

2.2.4. The fates of stellar ring 

The secular evolution of the gaseous disk leaves a stellar 
disk or a stellar ring for high and low viscosity, respectively. 
For a gaseous disk with high viscosity, the initial gaseous ring 
is moved to inner region through the star formation feedback, 
and forms a stellar disk with flexible width then, as shown 
in the £»-panel of Figure 2 and 3. It has been suggested by 
Noguchi (1999) and Elmegreen et al. (2008) that the bulges 
in local spiral galaxies form through the dynamical evolution 
of massive clumps of stars in high-redshift galaxies. The for- 
mation and growth of bulges may be related with fates of the 
stellar rings. When /? pea k < Rc, i-e. the stellar ring is at the 
outer edge of the bulges, the fates of the stellar ring might 
be determined by friction of the stellar ring with stars in the 
bulge since they are dense enough. In such a case, the stellar 
friction of the ring with stars in the bulge is driving the ring to 
immigrate into the bulge making growth of the bulges. This 
timescale is determined by 

^-0.3^) f^l.SxlO^iV^yr, (25) 

where t^ yn is the stellar dynamical timescale, /?! = /?/lkpc, 

V200 = Vrot /200km s" 1 , (To is the dispersion velocity (see eq. 
6 in Genzel et al. 2008) and V rot « ctq is used. We find that 
the timescale is significantly shorter than the typical episodic 
lifetime of star formation, namely, the ring immigrates into 



the bulge very fast. It happens for the case with high-viscosity. 
When the viscosity is low, the stellar ring will mixture with 
stars around its initial location. In such a case, the galactic 
stellar disk is growing. 

2.2.5. Dependence on £ 

For the linear star formation law, we have the critical value 
of q c in Paper I, which determines if there is a flat envelope of 
the Sg-evolution track. q c depends on the details of the index 
of star formation law, here we approximate q c ~ 1 .0 in light 
of the analytical solution. From equation (23), we have the 
critical value of £, 

r O 018 ^"1/3 -l/3 r -2/3 F _ 1/? -2/3 v 8/3 y (2-2 7 )/3 

U.U18a 02 q c C_ lQ h 5l K w V 3Q0 ^ 2 
Hc=\ (26) 
I 0.009 a-fq- c l E^V 3 2 m , 

for Case A and B, respectively. We refer the high-viscosity to 
£ > £ e whereas the low-viscosity to £ < £ e . We find that this 
simple estimation agrees with the numerical results as shown 
in Figure 2 and 3 for the cases of the high and low viscosity. 

It has been firmly indicated from observations that galax- 
ies are undergoing episodic star formation (e.g. Glazebrook 
et al. 1999). Kinetic feedback from SNexp in star formation 
depends on episodic active time (Afc). For an episodic activ- 
ity of star formation, the minimum stellar mass to produce an 
SNexp within Afc should be 

M c /M & > 7.0 Afo'j- 4 , (27) 

where we use the lifetime of hydrogen main sequence stars 

as f M s = 13 (M«/M Q )~ 2 ' 5 Gyr and Af .i = Af G /0.1Gyr. Oth- 
erwise stars less than this critical one will not play a role in 
the turbulent viscosity in star forming galaxies. Solution of 
the dynamical structure sensitively depends on the parameter 
£, which is 

£ = e/sN « w 1.0 x KT^wigf m^M^ , (28) 

where m .\ = M min /O.1M , m 7 = M C /7M Q and the Salpeter 
mass function with (3 = 2.35 is used. Here we show how the 
viscosity (£) depends on the IMF. Though £ oc Afg j 4 for 7M Q 
stars, we show £ oc At^ l for more massive stars in §5. Clearly, 
the kinetic feedback strongly depends on the IMF for high-z 
galaxies. 

We would like to point out that /sn can be in a quite large 
range in light of the IMF. For a top-heavy IMF, M m i n w 
0.5(1 +z) 2 M , where z is redshifts (Dave 2008), we have 
M min = 8M Q for stars at z = 3 and / SN w 4.6 x 10~ 2 Mq. On 
the other hand, for a flat IMF, (3 « 1.0 (Kroupa 2001), we 
have fsn = 0.125Mq. The /sn range is then expected to 
0.01-0.1Mq. For a given £, the larger /sn, the smaller e. 
The exact values of e and /sn are not important, but the £ with 
the similar range of /sn- 

2.2.6. Non-linear effects of star formation rates 

Compared with results in Wang et al. (2010), the non-linear 
star formation law used in the present paper generally leads to 
a dependence of the viscosity on gas density (see equations 
17 and 21). This non-linear effect results in a dependence of 
the structure and evolution of the gaseous disk on the initial 
conditions, especially on the initial surface density of the ring. 
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FIG. 5. — Left: Evolution of accretion rate of gas for Case A, here £ = 0.03Mg , corresponding to fig 2a. Lines with different colors (from red, green, blue, 
to pink) refer to different times from 10 6 yr to 10 9 yr with an interval A log? = 0.25. Right: the timescale of delivering gas with the star formation rates, the solid 
and dashed lines are for Case A and B, respectively. We find f a( jv ~ 10 s (SFR/ 10 2 MQyr~' )~°' 9 yr. 



For a linear star formation law, we find that the critical value 
of £ is £ c ~ 0.5C_8//?io significantly larger than £ c (equation 
26), where C-& = C*/10~ 8 and Rio = 7?o/10kpc is the initial ra- 
dius of the ring in Paper I, and is independent of the initial gas 
mass and density. For the KS law, 7 = 1 .4 leads to the depen- 
dence of q(R,t) on the surface density as shown by equation 
(23). Calculations are made for several cases with different 
initial surface density. We find that the gas surface density 
and stellar distribution for a non-linear star formation law are 
generally steeper than that of linear one, and the larger 7, the 
steeper E g and stellar distribution are. We do not provide the 
relevant figures here for saving pages. 

In a brief summary of this section, we obtain the numerical 
solutions of the gaseous disk and their properties. We draw a 
conclusion that a gaseous ring is going to form a stellar disk 
far away from it initial location through turbulent viscosity 
driven by SNexp. 

2.3. Delivering gas to galactic centers 

With the structure of the gaseous disks, we obtained in- 
flow rates for the issue as to feed SMBHs located at galac- 
tic centers. Figure 5 shows the inflow rates as a function of 
radius and time. We find that the resultant inflow rates are 
at a level of a few M Q yr _1 . A vanishing boundary condition 
for numerical convenience is used in this model at inner edge 
of the gaseous disk, the present model does not produce the 
true "central region". However, the model is enable to pro- 
vide the timescale of delivering gas and its dependence on the 
SFR. From §2.1.1, §2.1.2 and §2.2.1, we have the timescale 
favd = R/Vr, 



^adv " 



l/UR%T l) = 01 ^hEsKloVioo^Rll Gyr, 
=0.1 a^\ G EslC:{ Vf 00 ^ Gyr, 

(29) 



where £0.03 = £/0.03Mq for Case A and B. Here we only 
consider the high viscosity case, which is able to deliver gas 
into the center regions. Additionally, there are many evidence 



show that the IMF is top heavy in the circumnuclear star form- 
ing region that means high viscosity by SNexp in hundreds 
parsecs scale (see equation 28). The simple estimation of the 
timescale shows f ac jv oc £g 1,2 and f a d v oc £g 0,4 , where 7 = 1.4 
is used, for Case A and B, respectively, indicating that the 
timescale decreases with star formation rates. For a more 
clear impression, we use an approximation SFR ~ ttR 2 ^ as 
a global SF rates to show relation of f a dv versus SFR in the 
right of Figure 5. For typical star forming galaxies with star 
formation rates are of a few 10 2 M Q yr _1 , we find that the ad- 
vection timescale is less than 10 8 yr, namely, SMBH activity 
appears when starbursts precede f a d v < 10 8 yr for Case A. Case 
B shows much flat dependence of the time lag on the SFR. It 
should be noted that this difference is caused by the different 
dependence on S g . 

Transportation of gas in the ~ lpc scale (inside the dusty 
torus) in AGNs has been discussed by Wang et al. (2010). 
They show that the SNexp-driven inflows can form the 
Shakura-Sunyaev accretion disk around the SMBH with typ- 
ical sub-Eddington rates, naturally explain the well-known 
metallicity-luminosity relation. It is therefore that the deliv- 
ering timescale of gas is just the time lag of AGNs after star- 
bursts. We have to stress that the effects of thermal heating 
due to SNexp have not been included as feedback to the pro- 
cess of the star formation in the current model, which surely 
lower SF rates in the gaseous disk. This should be investigated 
by including of the thermal structure of the disk to improve the 
current model further. 



2.4. Multiple rings due to successive minor mergers 

The scheme described above is for a single episode of ac- 
tivity, however, recent observational results suggest that rapid 
but more continuous gas accretion via "cold flows" (Keres et 
al. 2005; Ocvirk et al. 2008; Brooks et al. 2009; Dekel et al. 
2009)) and/or minor mergers likely are playing an important 
role in driving star formation and mass growth of the massive 
star-forming galaxies at z > 1 (e.g. Noeske et al. 2007a,b; 
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Elbaz et al. 2007; Daddi et al. 2007£| . More interestingly, 
current observations of high-z galaxies show there is star for- 
mation inside of bulge with potential disk geometry (< a few 
kpc) while a star forming ring in the outer regions (Genzel 
et al. 2008). Successive minor mergers indicated by observa- 
tions are happening more frequently at high-z galaxies (Stew- 
art et al. 2009; Genzel et al. 2008; Forster Schreiber et al. 
2009). When the minor merger timescale (the inverse of mi- 
nor merger frequency) is longer than the advection timescale 
of gaseous disks, the multiple rings are evolving indepen- 
dently. In such a case, we only need include multiple minor 
mergers by an interval At for high-z galaxies. We do not 
include the interaction among the multiple rings. 

A bulge+ring structure has been often found in SINS sam- 
ple (Genzel et al. 2008). It may be due to multiple mergers 
in term of the merger rates with an interval of 10 8 yr or so 
(Lacey & Cole 1993). The intensive star formation taking 
place in the bulge is from the gas transported by the turbulent 
viscosity driven by SNexp and the ring or thick disk with star 
formation is due to a recent merger. For two successive mi- 
nor mergers, the first took some gas to the galaxy and made a 
gaseous ring, then it evolves and is transferred into circumnu- 
clear region in about 10 7_8 yr. The second minor merger will 
occur and create another ring structure. As a consequence of 
multiple minor mergers, the galaxy appears as ring+bulge. 



Molecular gas 
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FIG. 6. — The scenario of photionization of clumps hosting a star cluster. 
The geometry of the clump is enclosed, but the actual should be open like 
in M17 (their Figure 2 in Pellegrini et al. 2007) in order to transport kinetic 
energy of the SNexp into turbulence. The pressure balance is given by P lur = 
Pmd +^he an d Phg = ?^rad. where ? is around unity at the illuminated face (see 
text for a detail). 

For different minor mergers, the initial locations and the 
gas supply from the donor galaxies is different, leading to 
different behaviors of evolution and structure of the gaseous 
disks. As illustrations of the undergoing physical process in 
high-z galaxies, we try to apply the present model to two rep- 
resentative galaxies (BX 389 and BX 482) to show histories 
of evolution. We find that multiple mergers can set up the 
bulge+ring structure (see Figure 12 and Figure 14). We would 
like to stress that the multiple minor mergers are not arbitrar- 
ily adopted. On the contrary, the multiple model has been im- 

7 We would like to point out that the current model does not rely on sources 
of the gas supply whatever it is by minor mergers or cold flows. The present 
model does not predict a relation between merger rates and V r ot/Vtur. 



posed by more strong constraints: 1) the distance of the two 
successive rings, which is driven by the dynamical evolution; 
2) the interval of the two successive minor mergers, testing the 
merger rates; 3) the relative strength of star formation rates. 

3. PHOTOIONIZATION OF GASEOUS DISKS 

With the given structures of the evolving gaseous disk, we 
are able to, in principle, calculate the emission from the disk 
through photoionization by formed stars. For this goal, we 
have to specify following parameters: 1) ionizing continuum 
and luminosity; 2) the distance of ionizing source to the il- 
luminated face and geometry of photoionization; 3) density 
of the ionized gas and 4) chemical composition. Star forma- 
tion happens and forms a cluster mostly embedded in clumps. 
Each cluster containing a dozen to many millions of stars 
(Weidner & Kroupa 2005), which are surrounded by the lay- 
ers of the H II region, the photon-dominated region (PDR) 
and the molecular cloud. Generally, clusters are the energy 
source of ionizing their surrounding gas (e.g. Pellegrini et al. 
2007). Photoionization of clouds depends on several parame- 
ters, such as, the ionization parameters, geometry, density of 
the clouds and temperature. Following Lehnert et al. (2009), 
we use the ISM metallicity for all the radii of the gaseous 
disks. Figure 6 illustrates the geometry of clusters and physi- 
cal conditions of photoionization in one clump. We assume a 
closed sphere ionized by the center star cluster. 

We use the CLOUDY^ to calculate the emission from the 
gaseous disk. The photoionization is mainly determined by 
the ionization parameter defined by 



U(R,t) = 



$h(*,0 
n^(R,t)c' 



(30) 



where $h (cm s 1 ) is the surface flux of ionizing photons, «h 
is the total number density of hydrogen and c is the light speed 
(Osterbrock & Ferland 2006). Dopita et al. (2000) and Kew- 
ley et al. (2001) carry out detailed calculations of photoion- 
ization of star forming galaxies, but they took the ionization 
parameter Wasa free parameter. Unlike these papers, we cal- 
culate the ionization parameter in term of the given structure 
of the gaseous disk for photoionization. 

3.1. Ionizing spectrum 

To calculate the emission from the disk, we divide it into 50 
circular belts from 1 to 15 kpc. The surface gas density and 
star formation rates are given for each belt by the dynamical 
equations numerically. The status of the photoionized gas in 
each belt is determined by the radiation fields from the cluster. 
The ionizing spectrum, in principle, is determined by the IMF, 
the upper mass cutoff, cluster age and star formation history, 
etc. In this paper, we use the GALAXEV (Bruzual & Char- 
lot 2003; hereafter BC03) to calculate the ionizing spectra of 
young stars in star formation regions as the incident spectra. 

Given the initial mass function, metallicities and star for- 
mation history, the ionizing photons from a stellar population 
is given by 



1 



Ly(R, t) 

hv 



(31) 



where L v (R,t) is the total spectral luminosity at time t and 
radius R, and /?; on is the radius of the illuminated face of the 

8 http://ferland.org/cloudy/ 
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FIG. 7. — The upper panels show star formation histories at different radii for Case A. The lower panels are the ionizing spectra, which correspond to the star 
formation rate history in the upper panel. The initial radius of the gaseous ring is at R = lOkpc. The gaseous disk structure used corresponds to Figure 2a. 



TABLE 2 

Values of parameters in calculations 





Case A 






CaseB 




£(xlu 2 ) 


^0 


lo 


?(xl0 2 ) 


^0 


10 


3.0 


171.1 


0.70 


1.0 


133.56 


0.90 


2.0 


50.69 


2.36 


0.5 


66.78 


1.79 


1.0 


6.34 


18.87 


0.05 


6.68 


17.9 



We obtain dynamical structure of gaseous disk through con- 
tinue medium equations, but we have to connect the aver- 
aged structure with clump structures. Following Vollmer & 
Beckert (2002), we assume the turbulence with a spectrum of 
E(k) oc k~ 2 , where k is the wave number of the turbulence. The 
clumps have a characterized size comparable with the turbu- 
lent length (l c ~ /?sn)- We have the characteristic time scale 
that turbulence crosses the clumps 



ionized shell. We should note that L v (R,t) relates with the 
stellar evolution and is given by GALAXEV code. In the cal- 
culations, we take the minimum star mass M m i„ = O.IMq and 
the maximum stars M max = 100M Q . We use basic SSP models 
m62 (solar metallicity, Padova 1994 track) and Chabrier IMF 
for the ionizing source. 

With the typical values of the parameters given in Table 
2, we calculate the ionizing spectra from the star formation 
at different radii and times, which are shown in Figure 7 
and 8 for Case A and B at R = 2,4,6,8, 10, 12kpc and t = 
10 6 ,10 6 ' 5 ,10 7 ,10 7 ' 5 ,10 8 ,10 8 ' 5 ,10 9 years, respectively. The 
parameter a is fixed for each model in this section. The ion- 
izing fluxes just follow the star formation rates. The Chabrier 
IMF is used, the spectra thus keep a constant shape unless at 
the early stage of star formation and at the stage with a dra- 
matic decrease. The differences of the ionizing spectra for 
Case A and B result from the dynamical structure and the star 
formation history. 

3.2. ionization of clumps 



ftlir — 



(32) 



where l c is the characterized scale of th e clu mps, is equal to 
the gravitational free fall time = 1 / y/Gp c , where p c is the 
clump density and G is the gravitational constant. The char- 
acteristic scale of clump is approximated by l c = aH (Pringle 

1981). We have p c = n c m p = G~ l (V tm /aH) 2 



1.5 x 10 3 a ££tCZ 2 l0 E^R-*V^l- 2 ~<cm-\ 



0.2 xlO 3 a^^cm" 3 , 



(33) 



for Case A and B, respectively. Equation (33) shows n c oc 
Sg° 8 for Case A, n c tends to very large for two side edges of 
the ring since E g tends very small there. It is then expected 
that the two edges of the n c could be unphysical. Clusters 
within H II usually have a typical size of 0.1 to 10 pc, which 
is much smaller than the height of the gaseous disk at radius 
R. We introduce the filling factor C, which is defined by C = 
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Case B: 




X (A) 



FIG. 8. — Star formation histories and ionizing spectra as Figure 7, but are for Case B. The gaseous disk structure used corresponds to Figure 3a. 



g/Pc 



c = 



3.8 x 10" 3 (4_ 2 R W V^ 2 , 



6.2 x 10 a a2 £_ 2 t 5l C_ 10 L 2 K w V : 



10 K 300 ' 



(34) 

for Case A and B, respectively, where p e = T, & /H is the aver- 
aged mass density of gas as in continue disk and p c is the mass 
density of clumps. We will see C ~ 10~ 3 (in §4) if the density 
of clumps is determined by observations of [S II] A67 1 6 , 673 1 
lines (Lehnert et al. 2009). We note that the averaged ratio 
of cloud size to distance among clouds is about C 1 / 3 ~ 0.1 
which is enough for the validity of the fluid approximation in 
equation (1). 

We further specify the geometry of the photoionized clump 
as illustrated by Figure 6. Stellar winds, radiation and super- 
nova produce a cavity of hot gas inside the clump by pushing 
the ionization front outward (e.g. Pellegrini et al. 2007;Oster- 
brock & Ferland 2006). We take the wind-nebular interacting 
surface as the illuminated surface with a radius Ri on . If the 
radiation pressure from stars (P ra d) is given, we have 



(— 



1/2 



(35) 



where L is the luminosity of the stars in the cluster. In the fol- 
lowing, we make an attempt to build up the physical connec- 
tion between the photoionization process and the dynamical 
structure of the gaseous disk. 
The turbulent pressure balances the total internal pressures 



of the hot gas thermal pressure and radiation field pressure 
from stars, namely, P tul = Phg + Aad, where Phg is the thermal 
pressure of the hot gas. In principle, we have to calculate 
the hot gas pressure in term of stellar evolution and winds, 
however, it is very complicated and beyond the scope of the 
present paper. For simplicity, we introduce a parameter q, 
-fta = C^rad. we have 



^rad — fPtm j 



(36) 



where / = 1/(1 +?), P^ is given by the structure of the gaseous 
disk. Though the exact value of q is not known, it should be 
around unity ? ~ 1, otherwise the photoionzation is very weak 
for > 1 or the cavity is undergoing formation through the 
strong radiation pressure from stars for q -C 1 . We then we 
have / ~ 0.5. The photoionization is linked with the dynami- 
cal structure of the gaseous disk in light of equation (36). We 
have the ionization parameter 



U- 



P, 



rati 



/- 



h{v)nnc h{v)nyi ' h{v) PcC 11 



(37) 



where Fi 0n = Li on /47r/?i on is the ionizing flux, L loa is the ion- 
izing luminosity, m p is the proton mass, h is the Planck con- 
stant, p c = «H'«p and (v) is the average frequency of ionizing 
photons. The approximation is only used to clarify the de- 
pendence of U on /, rather than used in calculations. fj on is 
obtained from the stellar spectra in calculations. We use a 
constant density of gas throughout the clump. From the illu- 
minated face, the ionization structure of the shell can be self- 
consistently determined by heating and cooling in CLOUDY 
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for the given U at the illuminated surface. 

'Case A 1 




R (kpc) 



R (kpc) 



FIG. 9. — The photoionization models for given gaseous disks. From the 
top to bottom panels are ionization parameter U (at the illuminated face), 
Ar'Ptur, column density Nn (= 0.5S g /m ; , for photoionzation model), number 
density of clouds n c , the filling factor C and the Ha brightness surface S\\ a 
of the gaseous disk. Lines with different colors correspond to the times in 
Figure 7 and 8. 

It should be noted that the present ionization parameter 
U does not depend on star numbers (AQ of the cluster if 
since c or / are independent of N*. The hot gas pressure 
reads Phg = "hg£7h g , where k is the Boltzmann constant, «hg 
and 7h g are the density and temperature of the hot gas, re- 
spectively. Stellar wind supplies the hot gas through heat- 
ing of shocks produced by the interaction between the winds 
and medium of the inner surface. The power of the stellar 
winds is given by L w = E w /Y w , where E w ~ Af*M w V 2 , M w is 
the mass of stellar wind, V w is the wind speed and f w is the 
timescale of the wind replenishment. Considering the bal- 
ance of the shock heating and free-free cooling, f w should 
be equal to the free-free cooling timescale (%). This yields 

f w = fff oc n^T 1 / 2 , and L w oc Af*M w V 2 n hg r hg 1//2 ', which is bal- 
anced by the free-free cooling oc njLTig. We then have the 

9 -3 II 

hot gas density «hg oc iV*M w V^7jj and its thermal pressure 
P hg oc n hg T hg oc N*M W V^T^ /2 oc N^L\ J V^;T^ /2 , where we use 

M w oc L\ 1 (Pellegrini et al. 2007) and L* is star luminosity. 
The radiation pressure reads P rac j oc Af*L* from N* stars. It is 
expected that the parameter c; = F\i g /Pmd oc L% 1 is independent 
of N*, but dependent of inside the clump. From the simple 



argument, U is independent of N* 
and the density of the clumps. 



but depends on the IMF 



3.3. Emission from the gaseous disks 

We have 50 histories of star formation rates for these 50 
belts. This can be justified by the condition 

?adv >min(?Ms,fsF), (38) 

where fMS is the lifetime of the hydrogen main sequence stars, 
and ?sf is the timescale of star formation. We find that 50 
belts are accurate enough for the current calculations. For the 
./'-belt at time t , the ionizing luminosities are calculated from 
the star formation rates as 



AMJRj,t) = 2irRjARt+(Rj,t), 



(39) 



where Rj is the radius of the j-belt with a width AR. We 
get the total spectroscopic luminosity of the each belt from 
BC03, i.e., AMJRj,t) -t AL*JRj,t). Given the ionizing lu- 
minosity, ionization parameter and density of the ionized gas, 
CLOUDY produces the line luminosity ALi(Rj,t). The sur- 
face brightness of emission lines is given by 



Si(R h t) = 



AL,(Rj,t) 
2nRjAR ' 



(40) 



This equation produces the relation between ratio of emission 
lines and Sua, which can be compared with observations to 
test the theoretical model. We list all the necessary equations 
for the calculations of photoionization of the gaseous disk. 
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FIG. 10. — The upper panel shows the relation of the line ratio of [N II]/Hq 
with the Ha brightness. The lower panel shows the relation of the turbulence 
velocity and the brightness. Loops with different colors are obtained along 
the radius of the gaseous disk at the times corresponding to that in Figure 7 
and 8. 
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FIG. 1 1 . — The BPT diagram for the dynamical structure given by Figure 2-3 and 7-9. Left: Lines with different colors are at the same time in Figure 7 and 
8. Loops with different colors are obtained along the radius of the gaseous disk at the times. Right: stars with different colors correspond to the line ratios 
integrating the whole galaxies show the evolutions with time. It is clear that the galaxies are evolving along the arrow marked in the BPT diagram. 



3.4. Numerical results 

Figure 9-11 show the photoionization models correspond- 
ing to the dynamical structures for Case A and B. The ionizing 
fluxes follow the star formation rates of the gaseous disk. Fig- 
ure 9 shows the parameters of photoionization model along 
the radius. The ionization parameter U displays a synchrotron 
evolution with disk's dynamics, but holds a constant over the 
gaseous disk. The ionizing photons are shown in Figure 7 
and 8. We find that the gaseous disks are quite clumpy at a 
level of C ~ 10~ 3 . Ha brightness is calculated from equation 
(40) according to the photoionization model. The Ha ring is 
evolving following the evolution of the dynamical structure 
and star formation rates. 

Figure 10 shows the relation between line ratio [N Il]/Ha 
and Ha brightness in the upper panel. This plot actually re- 
flects the relation of the dynamical structure and photoioniza- 
tion. Clearly the trend is consistent with the observations as 
shown in Figure 1 1 in Lehenert et al. (2009). Emission line 
brightness is proportional to density and ionization parameter, 
namely, Si oc nJJ while the line ratios depends on the U~ 1 ' 2 . 
It is thus expected that there is an anti -relation between Si and 
line ratios for given the density. Star formation rates roughly 
determine the ionization parameter as well as the gas pressure 
inside the clouds (equation 36). Temperature is thus high as a 
consequence of the photoionization shown by the £ _I P-panel 
in Figure 9. The model produces a trend of turbulent veloc- 
ity with Ha brightness as shown in the lower panel of Figure 
10, which is consistent with Figure 6 in Lehnert et al. (2009). 
Generally, the brighter Ha ring the more intensive star forma- 
tion, and hence stronger turbulence. 

We calculate emission line ratios for the BPT diagram 
(Baldwin, Phillips & Terlevich 198 if] to test if there is evo- 

9 There are extensive calculations of the BPT diagram by adjusting the 
ionization parameter U (e.g. Dopita et al. 2000, Kewley et al. 2001), but there 
are no calculations combined with the dynamical structure of star forming 
galaxies. According to the results from Dopita et al. (2000) and Kewley et 
al. (2001), the line ratio relations are mainly determined by the ionization 
parameter hi. 



lutionary trend of star forming galaxies. In principle, the 
line ratios depend on the ionization parameter and metallicity 
mainly. Evolution and gradient of metallicity are complicated, 
but we fix the metallicity for simplicity to show how the line 
ratios evolve with ongoing star formation. We find that the 
ratio of [O m]/H/3 is decreasing with [N Il]/Ha, [S Il]/Ha 
and [O l]/Ha in term of increases of U. Figure 11 shows 
the BPT diagram of the gaseous disks for a series of given 
times. Since U is a function of radii as shown by Figure 9, the 
line ratios show open loops in the BPT diagram. Consider- 
ing the open loops are significantly overlapped, we integrate 
the whole galaxies to get the global ratios for the evolutionary 
trend of the galaxies. Obviously, the galaxies are evolving as 
shown in right part of Figure 1 1 . 

It is worth examining the evolutionary trend presented here, 
but it needs a large sample which allows us to estimate U 
through star formation rates. Starting from SDSS data, we 
group the sample in light of U and replot the BPT diagram to 
test the evolutionary trend. The results will be carried out in 
a separate paper. We would like to point out that we fix the 
metallicity for all radius and duration of evolution of galaxies. 
The real evolution of galaxies should be more complicated 
than that of pure star formation. It would be very motivated if 
we incorporate metallicity evolution with dynamics. 

3.5. Uncertainties of the present model 

The main uncertainty of the photoionization is the param- 
eter <, (or /), which appears in the ionization parameter U 
in equation (37). The ionization structure should be actu- 
ally self-consistently determined by the dynamics controlled 
by the stars and stellar evolution inside the clumps, however, 
this is extremely complicated and beyond the scope of the 
present paper. We invoke the observational constraints from 
Ml 7 cluster (Pellegrini et al. 2007) to simplify the complex 
problem. We show <; oc 7 , implying that U depends on IMF 
through both <^ and (v) . Though the uncertainties of this point, 
the main characters of the photionized structure of the gaseous 
disk remain. Future work on the structure and photoionization 
is worth exploring. Finally, we do not include the role of the 
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FIG. 12. — Evolution of for BX482 calculated from both Case A (right) and B (left) models. The parameters used are listed in Table 3. The thin black and blue 
lines represent evolution of the first and the second merged rings, respectively. The red lines are the current status of the rings at the observational times. The 
second ring evolves 3Myr since the minor merger. The interval of the two successive rings is ~ 10 8 yr. The profile of turbulent velocities agrees with the observed 
rings (Figure 3 in Genzel et al. 2008). 



mechanism energy of SNexp in line emission through shock 
heating as argued by Kewley et al. (2001). The role of shock 
heating in line emission mainly depends on the efficiency (£). 
We will consider the role of SNexp heating for line emission 
in the future paper. 

4. APPLICATIONS: BX 482 AND BX 389 

A few dozen of high-z galaxies are well observed currently 
in the SINS sample. They are often less evolved. As applica- 
tions of the present model, we choose two galaxies from SINS 
sample observed through the integral-field spectroscopy. BX 
482 and BX 389 are selected for applications, whose dynam- 
ical structure and emission lines are presented by Genzel et 
al. (2008) and Lehnert et al. (2009), respectively. BX 482 is 
rotating star-forming outer ring/disk dominated galaxy while 
BX 389 is central bulge/inner disk dominated (Genzel et al. 
2008). On the other hand, they are located in the star form- 
ing area in the BPT diagram (Figure 9 and 10 in Lehnert et 
al. 2009). This sets up additional constraints on the details of 
both dynamical structure and emission of the two galaxies. 

Generally, turbulence velocity, Ha surface brightness and 
[S ll]A6716/[S II] A673 1 can be measured from observations. 
For the free parameters listed in Table 1, we first estimate 



the orders of each parameters used in the model. Gas sur- 
face density is mainly determined by the Ha surface bright- 
ness, we have the estimation of E e . The parameter £ can 
be estimated from equation (15) or (19) in term of V tur for 
Case A and B, respectively. Given the estimation of £ and 
n c , a can be evaluated by equation (33). The timescale of 
dynamical evolution due to SNexp-driven viscosity is then 
given by f a av = R 1 jv. Though some parameters degenerate, 
we try to estimate them separately. The averaged mass den- 
sity is roughly n c ~ 1200cm" 3 estimated from the line ra- 
tio of [S II]A6716/[S IIJA6731 values (Lehnert et al. 2009). 
Given the dynamical structure, we calculate the photoioniza- 
tion model. 

4.1. BX482 

BX 482 is a ring-like star forming galaxy with small bulge 
(~2x 1O 1O M ) shown by Ha image, which has a ring radius 
of Rung = 7.0 ± 0.8kpc with a width AR ~ 1 kpc and a turbu- 
lence velocity V tul - = 100km s" 1 (Genzel et al. 2008). The 
galaxy is undergoing intensive star formation with a ratj^l 

10 The uncertainties of the star formation rates is larger than 50% in SINS 
sample (see Table 2 in Genzel et al. 2008). Here we take the uncertainties of 
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of ~ 140 ± 70M Q yr _1 . It shows two Gaussian profiles for 
the bulge and ring. The rotation curve measured by Ha line 
show a similar shape of rotation curve with equation (11). We 
set R c = 3.5kpc and V c = 300km s -1 for the rotation curve in 
this object. The star formation rate in the small bulge is of 
7^-buige ~ 35 ± 17M yr -1 (potentially with a disk geometry) 
while the star formation rate is of 7^disk ~ 105 ± 5 1M yr" 1 in 
the ring (estimated from Figure 3 in Genzel et al. 2008). This 
indicates that star formation is still ongoing in both the bulge 
and the disk of the galaxy yet. Clearly, it is necessary to have 
two successive minor mergers to explain the evolution of BX 
482. The initial parameters of the two rings are given in Table 
3. 



TABLE 3 

Values of parameters in applications 
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FIG. 13. — Photoionization model for BX 482 for the given structure in Fig- 
ure 12. The black and blue lines represent the current status of the evolving 
rings. The dotted and solid lines are for / = 0.4,0.6, respectively. We point 
out that the density n c oc E^ s 8 as shown by eq. (33) for Case A whereas the 
two lines overlap together for Case B. 

Figure 12 shows the dynamical structure of the models of 
Case A and B for BX 482. Generally, the models can explain 
the observed dynamical structures. £ for the second ring is 
smaller than that for the first in term of the smaller turbulence 
velocity. We integrate over the two rings and have the star 
formation rates. For Case A, we have 7£.buige ~ 33M Q yr _I in 
the bulge and 7?-disk ~ 92M Q yr _1 in the galactic disk. This 

50% for both BX 389 and BX 482. 
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agrees with the observed star formation rates. For Case B, we 
have T^buige ~ 17M©yr _1 and TZ^k ~ 90M Q yr _1 , which are 
still consistent with the observed within the uncertainties. The 
stellar surface £*, composed of the most old stars, is obtained 
from the integration of the star formation rates. The height of 
the gaseous disk is of 2kpc in the late phase of the evolution, 
which is rough in agreement with the scaled height observed 
(~ 1.6kpc). The profile of turbulence velocity of Case A is 
similar to the observed (compared with Figure 3 in Genzel et 
al. 2008) and better than that of Case B (which is too flat). 
Clearly the clouds of the gaseous disk are self-gravity domi- 
nated in light of Q < 1 driving star formation. The first ring 
has been immigrated into the bulge, which is still growing. 
The fate of the second ring is uncertain. It may stay at the 
current radius form a stellar disk for the given £ in Table 3. 
If £ is time-dependent and increases with time (discussed in 
§5), the ring is estimated that the second ring will immigrate 
into the bugle in about ~ 10 8 yrs as shown by the first and the 
third panel of Figure 12. 

Figure 13 shows photoionization models and emission from 
the star forming disks at the evolutionary epoch of 10 8 yrs. 
The two gaseous rings have very different U since they are at 
different evolutionary epoch as well as other physical param- 
eters. We find the gaseous disk is very clumpy as shown by C, 
which is of 10~ 2 ~ 10~ 3 . Equation (40) can produce the Ha 
emission which allows us to compare with Ha image. The last 
panel shows the surface brightness of Ha emission. We find 
that the results are in agreement with observed as shown in the 
distribution of the surface brightness (see Figure 3 in Genzel 
et al. 2008). We also plot the line ratios in the BPT diagram 
in Figure 16, which is consistent with the observations. 



4.2. BX389 

As another illustration, we apply the present model to BX 
389, which is a bulge-dominated star forming galaxy with 
a weak ring or disk at a scale of R^ng = 4.4 ± 0.5kpc and 
a turbulence velocity V tur = 180km s" 1 . This object is very 
different from BX 482. The total star formation rates are 
15O±75Af yr _1 , of which - lOO±5OM yr _1 in the bulge 
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FIG. 14. — Dynamical structure and evolution of BX 389, which show bulge-dominated structure. The lines have similar meanings but for different initial 
conditions as in Figure 12. The interval of the two rings are l() 8 yr. The turbulent velocity profile can be compared with Figure 7 in Genzel et al. (2008). 



and ~ 50±25M Q yr _1 in the ring. We set R c = 1.8kpc and 
V c = 315km/s to describe the BX 482's rotation curve (see 
Fig. 1). Other parameters are listed in Table 3. Two rings 
are assumed in modeling the dynamics and photoionization 
images. Since the star formation rates in the bulge dominate 
over in the weak ring or the disk (Genzel et al. 2008) and tur- 
bulence velocity in BX 389 is faster than that in BX 482, this 
implies that the dynamical evolution of the gas ring should 
be quite fast, namely, the parameter £ should be relatively 
higher than that in BX 482. This can be justified by equa- 
tion (15) or (19). We set the initial mass of the first ring is 
Mo,i = 5.6 x 10 10 M Q at radius ^o.i = 5.0kpc, which is mainly 
determined by the present star formation rate in the dynami- 
cal time. The black lines describe the dynamic evolution of 
the first ring. We find that the ring fast evolved inside the re- 
gion of the bulge with a timescale of 10 8 yr since £ is higher 
than that in BX 482. The second ring is injected after 0. lGyr 
of evolution of the first ring. The initial mass is determined 
by the surface density of Ha, but the dynamical age is un- 
known. The main difference from the first ring is the ini- 
tial mass, and other parameters are quite similar. Figure 14 
shows the dynamical structure and evolution of the two rings. 
The star formation rates in the bulge and disk are estimated 
as 7tbuige ~ 75M©yr _I and 7?.disk ~ 27M Q yr _1 for Case A, 



respectively, and T^buige ~ 70M Q yr _1 and 7^-disk ~ 27M©yr _1 
for Case B, respectively. The ratio of star formation rates is 
roughly one third in weak disk to the bulge, consistent with 
the observed. We note that the first merger ratio is about 1 :2.5, 
which is slightly larger than the definition of a minor merger. 
This could be called a small major merger, which still validate 
the present model approximately. 

Photoionization model is shown by Figure 15, which corre- 
sponds to the dynamical structure in Figure 14 for both Case 
A and B. The line ratios from the photoionization model are 
plotted in Figure 16 as the BPT diagram. The current model 
can produce the main features of dynamical structure of the 
two rings. The uncertainties of the model result in the dynam- 
ical age of the thick disk or weak ring, but the future track 
of the ring is predicted by the theoretical model. The spec- 
troscopy for the weak ring is not separated from the total ob- 
servationally. 

In summary, we show that emission and dynamical struc- 
ture are simultaneously reproduced by the present model. The 
SNexp, as a key ingredient in the model, is driving the evolu- 
tion of galaxies. It should be noted that the merged masses in 
both BX 389 and BX 482 are at a level of 10 I0 M Q . Such large 
mass accretions (> 10 I0 M Q ) must be accreted over roughly 
~ 10 8 years, which is consistent with the gas accretion rates 
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onto galaxies at these masses. If they are accreted on time 
scales shorter than this, they must be due to mergers, not cold 
gas accretion. The two galaxies favor mergers. 
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FIG. 15. — Photoionization model for BX 389 for the given structure in Fig- 
ure 13. The black and blue lines represent the current status of the evolving 
rings. The dotted and solid lines are for / = 0.4,0.6, respectively. Similar 
properties of n c can be found in the caption of Fig. 13. 



5. DISCUSSIONS 

We have presented detailed calculations of dynamical struc- 
ture and emission of the gaseous disk. Both rough estimates 
and numerical results of gaseous disk structures show that the 
parameter £ has strong influence on the structure and evolu- 
tion. In this section, we briefly discuss several potential fac- 
tors influencing £, such as, the roles of galactic winds, stellar 
evolution and initial mass function. The delayed SMBH ac- 
tivities and cosmological evolution as further implications of 
the present model are also discussed briefly. 

5.1. Influence of galactic winds 

Galactic-scale winds generally exist (e.g. Chen et al. 2010) 
and are thought to play an important role in evolution of 
galaxies. The present star formation rate density is much 
higher than the critical threshold of 0.1M Q yr -1 kpc" 2 , above 
which is powerful enough to drive vigorous galactic-scale out- 
flows in local galaxies (Heckman 2003). The outflows are ex- 
pected to play an important role in rendering a starburst "self- 
regulation" (Dopita & Ryder 1994; Silk 1997) supported by 
evidence of blue wings in rest-frame optical emission lines 



(e.g. Nesvadba et al. 2007). Its influence should be incor- 
porated into equations (1-3), even in the energy equation (6). 
Furthermore, broad emission lines of z ~ 2 galaxies have been 
found, which show Ha FWHM of 1500km s" 1 (Shapiro et al. 
2009). Such a broad component is thought to origin from 
galactic winds. The influence of stellar winds on the gaseous 
disk results in the gas loss and hence affects on the star for- 
mation. Transportation of angular momentum of the gas by 
SNexp-driven turbulence is stressed in this paper, which is 
supported by triggering of supermassive black holes hosted 
at centers of galaxies (Chen et al. 2009, see also Hobbs et 
al. 2010). Obviously, galactic winds will significantly lead to 
losses of mass and angular momentum, and hence change the 
evolutionary track of the gaseous disks. Likely, the peak in 
Sg-track appears in a further outer region in the presence of 
strong winds, even disappears. 

The well-known metallicity-mass relation (e.g. Panter et 
al. 2008) could be a result of galactic winds since less mas- 
sive galaxies have shallower potential so to allow a stronger 
outflows. It is expected that outflows have influence on the 
photoionization through metallicity. 



5.2. IMF, stellar evolution and metallicity 

There is increasing evidence for the cosmological evolu- 
tion of IMF (e.g. Kroupa 1999; Dave 2008), implying that the 
£ may undergo cosmic evolution. Top-heavy IMF appears in 
high redshift and may lead to more violent evolution of high-z 
galaxies, and AGN activities. Cosmic evolution of galaxies 
should include the role of IMF evolution in the dynamic equa- 
tions. Additionally, the correlation between metallicity and 
galaxy mass indicates that more massive galaxies are metal- 
richer (Tremonti et al. 2004; Brooks et al. 2007; Finlator 
& Dave 2008), namely, SNexp rates may undergo at differ- 
ent levels in galaxy evolution in light of star formation effi- 
ciency depending on galaxy mass. Obviously, massive galax- 
ies already underwent many episodes of star formation in the 
past, resulting in metal-rich. The present paper focuses on the 
evolution of structure on a dynamical timescale of galaxies, 
though the IMF determines the metallicity in a complicated 
way (beyond the scope of the present paper), we do not in- 
clude the cosmic evolution of IMF and fix it for all radii in 
the present calculations. IMF determines the SNexp rates, 
which control the dynamical structure of galaxies, leading to 
the metallicity distribution in galaxies. As a result of pho- 
toionization, metallicity govern the location of the galaxies in 
the BPT diagram. Future work by including the evolution of 
metallicity will provide more strong constraints on theoretical 
model and unveil the histories of mergers. 

We simply assume £ as a constant with radius and time 
within a dynamical timescale of galaxies, but it should be a 
function of radius due to the stellar mass function and stellar 
evolution. For massive stars, the hydrogen burning lifetime 

is given by f M s = 4.5 (M*/40M Q ) °' 43 Myr. For early phase 
of galaxy formation, the minimum mass of stars should be 
greater than M c > we have £ oc Af^ after 

inserting M c into equation (28), where A?5Myr = Af/5Myr. 
These massive stars will dramatically change the very early 
phase of the gaseous ring. SNexp rates as a function of time 
change £ by a factor of a few, may significantly change the re- 
sults presented here. This is the reason why we use a smaller 
one £ for the second ring which is only of a few 10 6 yrs in both 
BX 482 and BX 389 (see Table 3). 
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FIG. 16. — BPT diagram, the green circles are observation datas for the two galaxies. We take the line ratios (green symbols) from Figure 9 and 10 in Lehnert 
et al. (2009), which were estimated in a square aperture of 0.8" X 0.8" centered on the peak of the Ha emission. The model points (blue and red symbols) are 
integrated in the same regions. 



5.3. Inside-out growth of galaxies 

It is well-known that inside-out growth of galaxy disks 
drives the gradients of averaged age of stars, colors and metal- 
licity (Bouwens et al. 1997), however, M33 has a truncated 
disk at R ~ 8kpc (Ferguson et al. 2007) and the outer disk 
show an age increase with radius (Barker et al. 2007a, 2007b; 
Williams et al. 2009), indicating that the real growth of galaxy 
disks could be more complicated than the canonic version of 
inside-out growth. The present model only deals with the sin- 
gle episode of a minor merger or cold flows, and predicts three 
different fates of the merged gas as shown in Figure 2a-2c and 
Figure 3a-3c. For a cosmological evolution of galaxies, their 
growth histories contain various populations of minor merg- 
ers or cold flows, and form composite gradients of age and 
metallicity. As interesting applications of the present model, 
histories of galaxy growths will be formulated from the pop- 
ulation synthesis of minor mergers or cold flows. 

We use a presumed rotation curve in the paper, namely, we 
separate the gaseous disk from the stellar components of host 
galaxies. However, the inside-out growth of galaxies is gov- 
erning the potential on cosmic timescale. We do not include 
the potentials of the gas itself and the newly formed stars. 
This is valid for a single episodic activity of star formation. 
For a cosmic evolution of entire galaxies, we should change 
the rotation curves with inclusion of the newly formed stars 
in the last episodic activity and then get the history of star for- 
mation in the entire galaxies step by step. We would like to 
point out the rotation curves (equation 11) can be replaced by 
other forms under different potentials of galaxies. The present 
numerical scheme can conveniently include this. Future appli- 
cations to the large SINS sample should modify the rotation 
curves for more complicated structures and emission. 

5.4. Coevolution of SMBHs and galaxies 



Delivering gas into the centers of galaxies by SNexp will 
lead to coevolution of supermassive black holes and galax- 
ies, which is established for local galaxies (Magorrian et al. 
1998; Richstone et al. 1998) and may undergo cosmic evolu- 
tion in high-z galaxies and quasars (e.g. Coppin et al. 2008). 
Star forming galaxies at high redshift must undergo violent 
evolution: gas assembly, star formation, feedback from star 
formation and mergers and black hole activity. As a conse- 
quence, galaxies are growing fast at this epoch and may suf- 
fer from AGN feedback, showing increase of AGN duty cy- 
cle with redshifts (Wang et al. 2006; 2008). Supermassive 
black holes are triggered by mergers (e.g. Roos 1985; Carl- 
berg 1990; Wyithe & Loeb 2002), especially minor mergers 
(Li et al. 2009) with a natural delay at a level of the transporta- 
tion timescale following a star burst (Chen et al. 2009; Paper 
I; Netzer 2009; Kawakatu & Wada 2008; Davies et al. 2009; 
Wild et al. 2010). A unified model for galactic evolution 
and AGNs is required to describe the coevolution, in which 
feedback from star formation and AGN are included. The 
//-equation statistically describes the coevolution of SMBHs 
and galaxies (Wang et al. 2009a), but it should incorporate dy- 
namics of galaxies and accretion through transportation of an- 
gular momentum driven by SNexp turbulence of young mas- 
sive stars. In such a scenario, a delayed inflow to the central 
SMBHs is triggered in this way (Di Matteo et al. 2005; Hop- 
kins et al. 2006; Chen et al. 2009; Netzer 2009; Kumar & 
Johnson 2010). We can in principle establish the dynamics of 
coevolution. 

6. CONCLUSIONS 

We have a detailed extension of Paper I to investigate the 
roles of SNexp incorporated into the dynamics, for struc- 
tures and emission of evolving gaseous disk with given po- 
tential of host galaxies. Dynamical properties of such a 
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gaseous disk is totally different from the classical gaseous 
disk with a constant viscosity, which simply diffuses into 
space. We show that SNexp can efficiently transport gas 
into more inner regions and form a stellar ring or thick stel- 
lar disk. SBHM activity is then triggered with a lag of 

— 10 8 (SFR/10 2 M s yr 1 ) ' yr relative to starbursts. A stel- 
lar ring or disk forms, or immigrates into bulges depending 
on the efficiency £. For the case with high viscosity, the stel- 
lar ring formed from the secular evolution of the gaseous disk 
will immigrate into the bulge with dynamical timescale due 
to stellar friction. Bulges and SMBHs are jointly growing for 
high viscosity case whereas only galactic disks are growing 
for low viscosity. The model is able to produce the observed 
structures of high-z galaxies, including height, velocity dis- 
persion and young stellar components. 

Emission from the photoionized gaseous disk has been cal- 
culated through the simple stellar synthesis. UV photons are 
efficiently photoionizing the gas surrounding the star clusters 
and showing emission lines, such as [O III], [N II], Ha and 
H/3. The model predicts relations of line ratios and Ha bright- 
ness, BPT diagram. We apply the models to SINS sample, 
yielding the main features of the high-z galaxies. The consis- 
tency of the theoretical models with the observational proper- 



ties sheds light of the roles of SNexp in galaxy evolution. 

Future work to improve the present model are extensively 
discussed, especially the stellar evolution, galactic winds, ini- 
tial mass function and the geometry of ionized clumps. We 
stress the factors of influencing the efficiency £. When it is a 
function of time and radius, the structure and evolution of the 
gaseous disk will be changed somehow. The present model 
for single episodic event can be extended for successive mul- 
tiple minor mergers or cold flows, whose initial conditions are 
different from each other. The composite models by includ- 
ing all episodes could explain the metallicity and age distri- 
butions in galaxies in the frame of the inside-out growth, or 
more complicated observational results. 

The authors are grateful to the referee for a very construc- 
tive and helpful report significantly improving the manuscript. 
The paper is dedicated to JMW's father, who unfortunately 
passed away in June during the completion of the revised ver- 
sion of the paper, and expresses the infinite grief. We ap- 
preciate the stimulating discussions among the members of 
the IHEP AGN group. The research is supported by NSFC- 
10733010 and 10821061, CAS -KJCX2- Y W-T03 , and 973 
project (2009CB824800). 
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APPENDIX 
NUMERICAL SCHEME 



The non-linear partial differential equations (18) and (22) can be solved numerically by predictor-corrector methods (Lu & 
Guan 2003). For the standard form of the nonlinear differential equation. 



' Q=k l (x,Y)^+k 2 (x,Y), (0<x<l, andf > 0), 



Y(x,0) = g(x), 



(0<xsC 1), 



(Al) 



[7(0,0 = ^,(0,^^ = ^00, (t>0), 



where /j,\(t) and fi 2 (t) are functions of boundary conditions. We divide the disk into J (j = 0, 1,...,J) rings with a width h = dx. 
For simplicity, we define two operators as 



5Yj = Y H -Y H , t ,Y j = l -( Yj+l2+ Y H ). 
For the j-th ring at time nr, we have the predictor 

1 , n+i / x Y- +2 -Y? 



h 2 



S 2 Yp =h ( Xj ,YJ) 1 — 1 +k 2 ( Xj ,Yj) , (j= 1), 



and the corrector 



- 2 5\YJ + Yf) = kl (xjjp) J-—L^{x h Y^) , (j= l,...,J- 1), 



2i ^(Yj + Y r) = k 1 [x J ,Y r ) ' " 
where r = dt is the interval of the time. The initial condition is given by 

Yf=g(jh), (; = l,2,...,y-l). 

The inner boundary condition is given by 

y "=Mi(«r). 

The outer boundary can be expressed in two ways for the predictor and corrector 

^" = M2(y), and ±.0 (Y^Y^) = „2 (y ) , 
respectively. Setting j = 7 in equations (A3) and (A4) at the outer boundary, we have 

1 



Yj 2 — Y? 



J -+k 2 (xj,Y?) 



1 / i \ V n+l — V / i \ 

^ 2 S 2 {Y? + Yn=h (xj,Yp>) Y -^^ + k 2 (xj^) ■ 



(A2) 

(A3) 

(A4) 

(A5) 
(A6) 

(A7) 

(A8) 
(A9) 



Combining with (A7), (A8) and (A9), we can obtain Yj and Yj-\ by eliminating Y J+ \ . We get 7- 2 from (A3) and put it in (A4) 

j 



for T" +1 , and so on through the alternative direction implicit (ADI) method. 



